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1. Introduction 

Being a multiparticle quantum system, a magnet can be considered on different 
levels of hierarchy: a quantum (microscopic) level and a classical (macroscopic) one. 
The quantum level is described by means of quantum electrodynamics, or by simpler 
models like the Hubbard model or the Heisenberg one. The most common model for 
the classical level is the mean field model. Dynamics of a mean field is described by the 
equations of Landau-Lifshitz type. 

Each model is suitable to describe certain phenomena. For example, the problems 
of formation of large-scale structures (domain walls, topological solitons, nonlinear 
magnetization waves and so on) are naturally investigated from a classical point of 
view. More tenuous problems, like renormalization of the order parameter according to 
a temperature or an effective interaction constant, require a quantum point of view pp. 

Here we start from the quantum level described by the Heisenberg model. 
In addition to the usual Heisenberg bilinear interaction —J(S n ,S m ), we consider the 
biquadratic one —K(S n , S m ) 2 . By many theoretical and experimental researches it was 
shown that the biquadratic interactions have significant effects on magnetic properties. 
For example, a new ordered state (a nematic state, with zero magnetization) occurs as 
a separate phase transition [2J. Note, that the biquadratic interaction can be taken into 
account only if a magnetic system has the spin s ^ 1. 

In this paper we propose a classical generalization of the isotropic Landau-Lifshitz 
equation corresponding to the Heisenberg model with biquadratic exchange interaction. 
A transition from the quantum level to the classical one is performed by the mean 
field approximation. The classical model can be interpreted as a Hamiltonian system 
on a coadjoint orbit of the unitary group SU(3). Therefore, we acquire an additional 
mathematical apparatus, which gives a significant advantage. 

The mean field approximation gives a qualitative analysis of ordered states [3Pj . but 
has no answer about their stability. Moreover, in this approximation the temperature 
dependencies of order parameters considerably differ from the observed dependencies. 
That proves a necessity to take into account fluctuations of the mean field. The proposed 
effective classical models describe large-scale (or slow) fluctuations of mean field. 
One can come to slow fluctuations by an averaging over high frequencies pp. However, 
remaining in the context of theory of magnetism, we choose the models associated with 
the equations of Landau-Lifshitz type. 

The paper is organized as follows. Section 2 is devoted to the quantum model 
based on the spin Hamiltonian with biquadratic exchange interactions. We consider the 
SU(3)-invariant case. In section 3 we construct two effective models that describe large- 
scale fluctuations of mean field (the field of magnetization and quadrupole moments). 
We obtain one of them by an averaging of the quantum Hamiltonian over coherent states. 
The other effective model is a result of an averaging over mixed states. These classical 
models appear to be Hamiltonian systems on coadjoint orbits of the group SU(3), that 
follows from SU(3)-invariance of the original quantum model. Each coadjoint orbit is 
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determined by constrains, which are observed quantities becoming rigid after averaging. 
In section 4 we summarize results and give some ideas how to extend the proposed 
scheme to magnetic systems with higher spins. 

2. Quantum model of magnetic system 

2.1. Description of the model 

The magnetic system in question is realized by a homogeneous lattice of atoms with 
the spin s ^ 1 at each site. The lattice can be one-, two-, or three-dimensional, and 
has the distance I between the nearest-neighbor sites. We assign three spin operators 
S%) to each site n; they obey the standard commutation relations: 

rqa q/Si _ i^aP-yq-yr 

where a, (3, 7 run over the set {1, 2, 3}, and 5 nm denotes the Kronecker symbol. 

We use the localized spin model for the magnetic system. In many cases this 
model adequately describes a magnetic system by the Heisenberg Hamiltonian, which 
includes only the bilinear exchange interaction. Nevertheless, there are a lot of magnets 
that require taking into account higher powers of exchange interaction. Our model is 
applicable to magnets with the spin s ^ 1. 

In the present paper we consider the Hamiltonian with biquadratic exchange and 
call it bilinear-biquadratic: 

H = - ^{ J (Sn, S n+5 ) + K(S n , S n+S ) 2 }, (1) 

n,<5 

where S n = is a vector of spin operators at site n, and S runs over the 

nearest-neighbor sites. This Hamiltonian was discussed, for example, in [2H5]. The 
constants J and K serve as exchange integrals. We suppose that J and K are positive. 
It means that we consider a ferromagnetic interaction in preference. 

The operators {5""} (here n is fixed) are defined over the (2s + l)-dimensional space 
of irreducible representation of the group SU(2). They generate an associative matrix 
algebra over this space. The complete matrix algebra can be represented as a direct 
sum of irreducible sets of tensor operators with respect to the action ad^ Q . In the case 
of s — 1, we have: Mat 3x3 ~ [9] = [1] + [3] + [5]. Evidently, the operators {5*°} form a 
basis in the 3-dimensional irreducible set. One can construct a basis in the 5-dimensional 
irreducible set from the tensor operators of weight 2. These are the quadrupole operators 
{Qn, Qn, Qn, QJ» ,2] , Qn'^} defined by the formulas: 

The spin and quadrupole operators are normalized by the following relation: 
Tr(P) 2 = §s(s + l)(2s + l). 
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As s — 1 we have Tr(P) 2 = 2 . The chosen normalization is matched to the relation 

(^) 2 + (^) 2 + (^) 2 = ^ + l)- 

Now, fix the canonical basis {|+1), |— 1), |0)} in the space of representation. Then 
one obtains the following matrix representation for the spin and quadrupole operators: 
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We denote all spin and quadrupole operators: {S^, S^, S^, Q}?, Q}?, Q"^, Q 
by {Pn}a=i- The operators {P^} obey the following commutation relations 

[P n i P m ] iC a b c P n 5 nrn , 

where C a bc are structure constants; nonzero components are 

Cl23 = C*145 = C*167 = C*264 = C*257 = C*356 = 1; 
Cl68 — C*528 — V^3, C437 = 2. 

The Hamiltonian ([1]) becomes bilinear in the terms of {P^Y- 

% = ~{J ~ \K) z2 ^n^n+S ~ \K 52 52 QnQn+5 



[2,2] ^[2,0] } 



n,<5 a 



-KN 

3 ' 



(2) 



where A denotes the total number of sites. Obviously, the Hamiltonian is SU(2)- 
invariant, and one can transform the operators {S%} and {Q^} by the formulas of 
adjoint representation 

US«U~ l = D a \U)S p n , D a ? e SO(3), 



UQIU- 1 = D ab (U)Ql D ab e SO(5), 



where D al3 (U) and D ab (U) are matrices of the real irreducible 3- and 5-dimensional 
representations of the group SU(2) respectively, and U— exp{^ Q ip a S®}, where {<p a } 
are group parameters. As K = J the SU(2)-symmetry is extended to the SU(3)-one, 
and the Hamiltonian (T5]) gets the form 

n,S a 
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2.2. Mean field approach and ordered states 

Instead of interactions between the spin and quadrupole operators {P®} according to the 
Hamiltonian (j2J), we consider effective interactions of the operators {P%} with a classical 
mean field. We suppose that components of the mean field at site n are proportional to 
averages (quasiaverages) of the quantum operators {-P^}. 

In the mean field approximation the Hamiltonian (JSJ) has the form 

^MF = —{J - \K)Z £ £ ^(^> - \KZ E Qn(Qn) ~ l KN *> ( 4 ) 

n a n a 

where z is a number of the nearest-neighbor sites. We have to give a warning about 
averages of {P"}. If one calculates the averages by means of the density matrix 
p(T) = exp{ — ^}, one obtains zeros. This follows from the SU(2)-symmetry of the 
Hamiltonian (j2J). Nonzero values of the averages appear if the symmetry is broken. 
Symmetry breaking can be stimulated by an external magnetic field that vanishes after 
specifying an order in the magnetic system. Such averages are called quasiaverages [BJ. 

Suppose that the magnetic system in question has nonzero quasiaverages {(P")}- 
They form a classical 8-component vector field {/j, a (x n )}^ =1 , which we call a mean field. 
Suppose that the mean field is constant over the whole magnetic system. This happens 
in the case of thermodynamic equilibrium and an infinite lattice. Then under an action 
of the group SU(2) the Hamiltonian (J3J) can be reduced to a diagonal form, namely: 

n n 

= {( J " 5*0^3 + iKQ^fis + Ik} , 

n 

where the components fi 3 = (S 3 ) and fi 8 = (Q' 2,0 ') do not depend on the spatial point x n . 
These components are suitable to be order parameters. Evidently, fi 3 describes a 
normalized magnetization (a ratio of z-projection of magnetic moment to a saturation 
magnetization), is similarly connected to a quadrupole moment. 

Now we briefly show that the proposed quantum model admits ordered states. In 
the mean field approximation a partition function is calculated by the formula 

/IMP 

Z(fi 3 ^ 8 ,T) = Tre~ , 

where /imf denotes the one-site Hamiltonian 

h UF = ~(J - \K)^S Z - \K^Q^ - IK. 

The mentioned mean field exists if self-consistent relations are held, in other words, if 
the system 
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has a solution. After calculation of the mean field averages one obtains the self-consistent 
relations in the form 



2 sinh 



,, _ kT 
^3 



2 cosh^-expj-^} 
V / 3exp{-i|M} + 2sinh^ 



Solutions of the system correspond to ordered states of the magnetic system in question. 

An evident solution is the paramagnetic state (/a 3 = 0, fi 8 = 0). All other solutions 
depend on a temperature T, and the exchange integrals J and K. Note, that we 
consider the ferromagnetic interaction in preference: J > 0. Nontrivial solutions appear 
at temperatures low than the critical one T crit = ^{J — \K). As K<0 there exists 
a ferromagnetic state with the values (/a 3 — 1, ^ — ^) at zero temperature, and a 
nematic state with the values (/x 3 = 0, /i 8 = ^=) at zero temperature. As i^>0 there 
exist four nontrivial solutions: two ferromagnetic states with the values (/i 3 = 1, = ^75) 
and (/i 3 = |, ^8 = ^Tg) a t zero temperature, and two nematic states with the values 
(/i 3 = 0, ^ = 5§) and (yU 3 = 0, /is = ^) at zero temperature. The same states are 
declared in [3111]. The states (/i 3 = l, ^8 = ^75) and (/x 3 = 0, A i 8 = ^7§) are stable. The 
problem of transient processes in the mean field approach is discussed, for example, 
in [1]. The analysis of solutions of the self-consistent relations proves that ordered 
states in the proposed model exist. 

In the sequel we deal with the case J = K, which corresponds to the boundary 
between the ferromagnetic and the nematic regions (see the phase diagram of the 
bilinear-biquadratic s = 1 model in [5]). In this case, the Hamiltonian fl2D and its mean 
field approximation are SU(3)-invariant. The latter gets the form 

n MF = -lJ Z J2J2 P n( P n) ~PNZ = -IJ Z J2J2 P >« ~ l JNZ - ( 5 ) 



2.3. Motion equations for large-scale fluctuations of mean field 

Return to the quantum SU(3)-invariant spin model with the Hamiltonian fl3]). The 
Heisenberg equation for an evolution of P° has the form 

ift^f = (6) 

We suppose that the magnetic system is ordered, then we take an average of equation (J6|) 
over the Heisenberg (time independent) coherent states 



|^H) = -^(ciH|l)+c_ 1 (n)|-l) + c H|0)), 



N 

|ci| 2 + |c_i| 2 + |c | 2 = 1. 
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Alternatively, one can take an average by means of the density matrix. In the both cases 
we neglect correlations between fluctuations of the quantum fields {Pn}a=i a ^ distinct 
sites, that is 



'i 

\ ~, / r>a\ I r>b ' 



(P a n P b m ) « (P a n )(P b m ) = fi a (x n )^ b (x m ). (7) 
An averaging of equation (jBj) results in the following equation for fx a (x n )\ 

fl = ^J^Cabc^biXn) \jl C>XX {X n ) + H C ,yy(X n )) , (8) 

which is a Hamiltonian one with respect to the Lie-Poisson bracket. 

In order to investigate large-scale fluctuations of the mean field {[i a (x n )}l =1 , we 
consider a continuum space instead of the discrete lattice. It can be achieved by 
the well-known limiting process. In the case of SU(2)-magnetic system (only bilinear 
intereations are taken into account), this limiting process underlies the macroscopic 
phenomenological theory of magnetism [7j. The limiting process replaces quantum 
operators by densities of their averages, which serve as dynamical variables. In our 
case, we deal with the densities M a of averages of the spin and quadrupole moments: 

M a (x) = Y ( (i a {x n )S{x,x n ), S(x,x n ) = \ J ^fjjfs 
^ y x n £U[x), 

where Vq denotes a physically infinitesimal region of the lattice, and U(x) is the 
infinitesimal neighborhood of x. The Lie-Poisson bracket for {M a (x)} is defined by 

{M a (x), M b (y)} = C abc M c (x) 5{x - y), 

where S(x) is the Dirac function. Since dimensionless quantities are more suitable, 
we introduce fi a (x) = VoM a (x) instead of M a (x). Then equation (jHJ) gets the form 

n ^- = {W. cS , fJ-a{X)} = VoCabcfl b (x)— , (9) 

at dfj, c 



where I is the lattice distance, and d is the lattice dimension. Note, that in the 
2-dimensional case we obtain a scale-invariant Hamiltonian. 

Evidently, (Q is a generalization of the well-known Landau-Lifshitz equation to the 
case of 8-component vector field {fi a }- In the same way one can obtain the standard 
Landau-Lifshitz equation, if considers a spin system with s = | over the 2-dimensional 
space of representation of SU(2). 

We rewrite OH]) i n the matrix form 

n tt =^rl&Afl, » = Y,^ pa - ( 10 ) 

a 

Here fi is a Hermitian 3x3 matrix, [-, ■] denotes the matrix commutator, A is the Laplas 
operator. Being SU(3)-invariant equation f JTUj) as well as © preserves the quantities 
ho = \ Tr fi 2 and fo = \ Tr p? , which we call invariants. They serve as constrains for the 
Hamiltonian system and define the manifold where the vector field {/i a } lives. At the 
same time, this manifold is an orbit of coadjoint representation of the group SU(3). 
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3. Classical Hamiltonian systems on coadjoint orbits of SU(3) 

In the 1-dimensional case the Hamiltonian system appears to be integrable, what is 
shown below by means of the orbital approach. 

3.1. Phase space for SU (3)- symmetric generalization of Landau-Lifshitz equation 

In this section we briefly construct the orbital interpretation of a finite-zone phase space 
for the SU(3)-symmetric generalization of the Landau-Lifshitz equation. 

Consider an algebra of polynomials in A with coefficients from the Lie algebra su(3). 
Denote by q + the algebra su(3) ® V(X), where P(A) is a ring of polynomials in A with 
the standard multiplication. Let A, B e g + have the form: 

N+1 N+1 

A(X) = ^ A ™' B ( x ) = E ^' Afc ' ^ n > B k e 5u(3). 



n=0 fc=0 



Then 



[A 1 B]=J2iA n ,B k }X n+k eQ + . (11) 

n,k 

The operation (TTTT) turns $j + into a graded Lie algebra. 

Let P a ' n = X n P a , where a runs from 1 to 8. The set {P a > n } serves as a basis in g + . 
Recall that [P a , P b ] = iC a b c P c ; the nonzero components C a b c have the following values: 

C123 = Cl45 = Cl67 = C*264 = C257 = C356 = 1, 

C*i68 — C528 — Vo) C437 = 2. 
Introduce a bilinear ad-invariant form on g + by 

(A,B) = ^resX- N - 2 TrA(X)B(X). (12) 

The basis {P a>n } is orthonormal with respect to the bilinear form. Let M. = jj^ be a 
dual space to the algebra g + with respect to ( TT2"]) . Orthonormality of {P a < n } implies 
that {p a > n } also form a basis in Ai. Consider the following elements of M.: 

N 8 

= E E ^ An/,a + (^ +1/?3 + ^ +lp8 ) AJV+1 - 

n=0 a=l 

The functions /2(A) form a closed ad-invariant subset of .M, we denote it by A4 N+1 . One 
can compute the coordinate /1™ of /2(A) by the formula 

^ = (A(A),P a '- n+w+1 ). 
Define a Lie-Poisson bracket in C(Ai N+1 ) as 

m,n a,b " 
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with the Poisson tensor field 

Introduce also two ad-invariant functions /2(A) and /3(A) by the formulas 
/ 2 (A) = ^Tr/i 2 (A) = ^^(A), 

a 

J 3 (A) = ^Tr/i 3 (A) = y /f |d o6c /i a (A)/i 6 (A)/i c (A), 

where 4&c = ^| Tr(P a P b P c + P b P a P c ), and // a (A) denotes the polynomial 
MA)=^ + ^A + ^A 2 + ---+/,r i A^ 1 . 
The invariant functions are also polynomials in A: 
I 2 (X) = h + h 1 X + --- + h 2N+2 X 2N+2 , 

hW = fo + flX + --- + f3N +3 X 3N+3 . 

It is easy to prove that the coefficients {ho, . . ., hjy+i, fo, . . ., /jv+i} are annihilators 
with respect to the bracket f lTHj) . We fix these coefficients and obtain the system of 
algebraic equations 

h n = const, f n = const, n — 0, . . . , N + 1, (14) 

which determines an embedding of an orbit O n+1 of dimension 6(iV+l) into A4 N+1 . The 
coefficients {h^ +2 , • • • , h 2 ^ +2 , In+2, ■ ■ ■ , ^3^+3} are pairwise commutative integrals 
of motion. We call them Hamiltonians. In the 1-dimensional case the number of 
Hamiltonians is sufficient for integrability of the Hamiltonian system on an orbit. 

Here we are interested in two Hamiltonians: h N+2 , h N+3 , and the corresponding 
Hamiltonian flows. The Hamiltonian /i/v+2 gives rise to the stationary flow 

{X, h N+2 } = 2C ab A n c + \ a = 1, . . . , 8. (15) 



dx 

The Hamiltonian h^ +3 gives rise to the evolutionary flow 

^ = {^,h N+3 } = 2C a Mtf +2 + dv n c +1 ), o = l,...,8. (16) 

Equations (ITS]) and f fTBj) are compatible, for the corresponding Hamiltonians 
commute: {/iat+2, /iat+3} = 0. Thus, (fT6|) describes an evolution on the trajectories of 
(IT5|) . that is the dynamical variables {/i™} in (ITBj) depend on x. From (|T5|) and (ITBj) we 
have: 

— - 2C a&c/ i 6/ i c - — . (17) 

The variables can be expressed in the terms of {fJ%} and then (|T7j) becomes 

a closed system of partial equations for In order to compute the variables {yU^} 

one has to solve the following degenerate system of equations of the stationary flow: 

d ^ = 2C ab M, o = l,... ,8. (18) 
It becomes possible if one restricts the system to the orbit Q N+1 C Ai N+1 . 
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3.2. Classification of orbits of SU(3) 

It is evident, that the orbit O n+1 defined by (|14p is a vector bundle over a coadjoint 
orbit of the group SU(3). That is why we need to classify orbits of SU(3). 

The group SU(3) is simple [8], hence its algebra jj~su(3) coincides with the dual 
space g*. Consequently, the coordinates {/^ a } in 0* can be regarded as coordinates in 
5u(3) as well as in su*(3). A generic element fi Gsu*(3) has the form 



// 



\ V2 



M3 + 
u 7 + i/x 4 

(/ii + i/t 6 + Ms + 1/^2} 



1M4 

1 



1 

v/2 



-M3 + 7^8 

(fit + i/t 6 - /ts - 1M2) 



1 

^2 
7^ 



l/*6 
__2_ 

V3 



M5 + i/*2) 



(19) 



/ 



Let I) be the maximal commutative subalgebra (also called the Cartan subalgebra) of g. 
The dual space f)* to the Cartan subalgebra f) coincides with I). 

By definition the set C^ in = {g _1 /ii n g, VgGSU(3)} is the coadjoint orbit of SU(3) 
through an initial point p,- m e 5U* (3) . All elements (?' G SU(3) such that g'~ x fi\ xi g' = /to 
form the stationary subgroup at /ti n . The orbit (9^ in is a homogeneous space, 
which is diffeomorphic to the coset space SU(3)/S^ in . There exist two types of orbits of 



SU(3) 
U(l)xU(l) 



of dimension 6, and the degenerate Odeg 



SU(3) 



SU(2)xU(l) 



SU(3): the generic O gcn 
of dimension 4. 

It is proven by R. Bott that each orbit of coadjoint action of a semisimple group G 
intersects f)* precisely in an orbit of the Weyl group W(G). 

The full Weyl group of SU(3) consists of six elements {e, <j\, cr 2 , o"i<7 2 , o" 2 o"i, cr 1^2^1)1 
where ai, 02 are reflections across the hyperp lanes orthogonal to the simple roots 
oti, a 2 (see figure 1). The open domain C—{ji&i)*, (ft, a) > 0, Wa G A + } is called a 
positive Weyl chamber. Here A + denotes the set of positive roots. We call the set 
To = {/t G f)*, (/t, a) = 0} a wall of the Weyl chamber. An orbit of the Weyl group W(G) 
is obtained by the action of W(G) on a point of C. 



Weyl chamber 



Figure 1. Root diagram of SU(3). 

Each orbit of the Weyl group W(G), and consequently, each coadjoint orbit of G 
intersects the positive Weyl chamber in an only point. That is why we can classify 
coadjoint orbits of G by points of the positive Weyl chamber. 

In the case of group SU(3), there exist two types of orbits of the Weyl group. 
A generic orbit contains six elements and passes through the interior of the positive 
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Weyl chamber. A degenerate orbit contains three elements and passes through a wall 
of the positive Weyl chamber. According to this, we call an orbit of SU(3) a generic 
one if /ti n lies in the interior of the positive Weyl chamber, and a degenerate one if jl- m 
belongs to a wall of the positive Weyl chamber. 

In our case, has the following diagonal form 

/2 in = diag(m + ± q, -m + ± q, q), 

where m and q denote initial values of the variables /13 and respectively, or boundary 
values (at zero temperature) of the corresponding components of the mean field. 
As m > 0, q > the coadjoint action of SU(3) gives a generic orbit. If m = or m = a/3 q, 
we obtain a degenerate orbit. In the sequel we consider degenerate orbits with m = 0. 

3.3. Hamiltonian equations on orbits of SU(3) 

Return to the system of equations ffl8l) . which is degenerate in Ai N+1 . However, it 
can be solved if one restricts the system to the orbit O n+1 C Ai N+1 . Each orbit is 
determined by the following equation [9]: 

Xmin(A) = 0, (20) 

where Xmin(A) * s the minimal characteristic polynomial in £iEO N+1 . Equation ( 1201) 
serves as a constrain for the system f llSjl . which has the form 

(9/0° 

-£=Ad fl0 »\ (21) 

Now we solve (I2T1) on orbits of the group SU(3). 

A degenerate orbit is determined by the equation 

^ + \pf V - IT = °> 
where ho = q 2 = const. Using this constrain, one obtains the following solution of ( 12~TT) : 
= 6^ Cabc^b^x + ^ A*a) where ^ /i° is an element of KerAd^o. The motion 
equation (TTTj) on the degenerate orbit has the form 

dfJ-a _ 8A r 1 SAM.. (nn\ 

Q t ~ 3h y ~ / abcH'bt J 'c,xx T Hq ^a,x, \^) 

where we write \L a instead of /x° and scale the flow parameter t by 16*4. The dimensional 
constant A provides a correspondence between (I2"2"j) as hi = and ffTUj) as d—1. That 
is ( 1221) describes large-scale fluctuations of the mean field {n a }- 
A generic orbit is determined by the characteristic equation 

/t 3 - h p, - §/„ = 0, 

where ho = m 2 + q 2 , and fo = -^(3m 2 q — q 3 ). On this orbit we obtain the following 
solution of (|2II»: 

^« = 8(hg-3/ 2 ) (hlCabc^b^x - 2V3 foCabcVb^x + h oC a bcVbVc,x) + 
2/q/i -3/tgfti Q , 3/ /n-2fto/i Q 
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where 77° = a/5 dabcfJ^^c- The motion equation (TIT]) on the generic orbit has the form 

~Qj~ = fig- ^3/^ (j^O^abc^bf^^xx V^3 foC"abc^bVc,xx 

— v3 foCabcTJbficxx + h Q C abcllbll^xx ) + (23) 



I 8.4 2/o/i-3/igfei 1 8.4 3/o/n-2/to/i 

"T 3 / 2 -fig /V* I" 3v /3 /o 2_ ft 3 '/a,x- 



As /ii = 0, fi = equation ( 123]) also describes large-scale fluctuations of the mean field. 
One can obtain f[2"3"]) from ([H]) by averaging with a more complicate correlation rule. 
Equations (1221) and (1231) imply the following Hamiltonians respectively: 

^ a=l 

%gen = ftff/g / y^faC/Va:) 2 + h (r] ajX ) 2 - 2\/3 foPa.xV^x J dx. 
a=l 

In addition to the 1-dimensional case one can consider the corresponding 2- or 
3-dimensional Hamiltonian systems with the effective Hamiltonians 

U cS = J J H(fi) d d x, (24) 

where // denotes {/x a }® =1 . The exchange integral J = A/h gives the Hamiltonian the 
required physical dimension. By H we denote the Hamiltonian density 

d 8 

#deg = 3^ J 2 ' OT 

fc=l o=l 

rf 8 

#gen = X|Z3^ ^2^2(ho(Va,x k ) 2 + h {r) a>Xh ) 2 - 2^3 f ^a,x k Va,x k j ■ 
fc=l a=l 

One can use these effective Hamiltonians for describing the magnetic system considered 
in section 2. Note, that Hdeg is the same as the Hamiltonian of ([9]). 

The proposed Hamiltonians describe large-scale (slow) fluctuations of the mean 
field /x. After averaging over high frequencies some observed quantities become rigid 
(or invariant); these quantities are /io = <Ww*&, an d /o= a/5/3 dabcHa^b^c- They serve 
as constrains for the Hamiltonian systems, and are equivalent to (|2"D1) . The constrains 
determine the orbit where the system has to be considered. 

In the case of SU(3)-invariant model we deal with the magnet whose ferromagnetic 
and nematic states are equiprobable. A generic orbit corresponds to a state with the 
ferromagnetic order at zero temperature, because of nonzero magnetization (m^O). 
A degenerate orbit (m = 0) corresponds to a state with the nematic order at zero 
temperature. So equations (122]) and ( 123]) describe fluctuations of the mean field /x 
near a nematic and a ferromagnetic ordered states respectively. 
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3.4- SU(3)-invariance of effective Hamiltonians 

As mentioned in Section 2, the quantum Hamiltonian ^ and the mean field Hamiltonian 
(jlj) are SU(3)-invariant as K=J. Here we show that the proposed classical effective 
Hamiltonians 023]) are also SU(3)-invariant. 

Recall, that the mean field {/x a } belongs to the real 8-dimensional space of coadjoint 
representation of SU(3). Hence, an action of SU(3) transforms {^a} by the formula 

f/, a = Dab^b, D ab e SO (8), 

D a b is a matrix of the real irreducible 8-dimensional representation of the group SU(3). 

Note, that the tensor d abc satisfies the relation d abc d qbc = S aq . The components {d abc } 
serve as Clebsch-Gordon coefficients for a decomposition of tensor square of the coadjoint 
representation into irreducible components. In this connection, we have the following 
relation, well-known in theory of representations, D^iD^j = d qbc d q iy c iD qq i . Then as a 
result of the action of SU(3) on {r] a } we get 

rj a = \fhd ahc b w \i hl b C( ,,\i d = b qq iT] q i . 

The action of SU(3) on the vector fields {[i a ,x} and {t] a ^ x } is the same. Therefore, the 
densities H gen and H^eg are SU(3)-invariant. 

Densities of the effective Hamiltonians can be expressed as 

jk ab ^ 

where g a b{^) serves as a metrics invariant under an action of the group that 
transforms /x, and Gjk(x) is a metrics in the a?-space. For the proposed effective 
Hamiltonians the cc-space is Euclidean: Gjk{x) = 5jk- The metrics in /j,-space is trivial: 
dabifJ') = 8 a b in the case of a degenerate orbit, and has a more complicate form: 

g a b{v) = h 3^ 3f 2 [hl^ab + 20h d cpa d cqb fi p [i q - 4 V / 15 f d abc fi c ^ 

in the case of a generic orbit. 

The density (125j) can be interpreted as a Lagrangian density of relativistic o- 
model; in this case Gjk is the metrics of the Minkowski space. After quantization one 
obtains a Hamiltonian system that describes slow fluctuations. Quick fluctuations can 
be taken into account by means of a renormalization group pQ . It makes the coefficients 
,3 \ f2 and 7T7~ dependent on parameters of the renormalization group, for example on 
a temperature. 

3. 5. Parametrization of orbits 

Remarkably, that the effective models are entirely defined by geometry of orbits. We 
will prove this statement, if perform a parametrization of orbits and express the effective 
Hamiltonians in terms of these parameters. 
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A generalized stereographic projection gives a suitable way of parametrization for 
coadjoint orbits of a semisimple Lie group [10]. In the case of group SU(3) we have 



where 



1 U)2+Ul 3 +Ul2+Ul 3 C — 1 (l-Wi)(w3-WiW2) + (l-Wi)(w3-WlW2) 

^2 l+\w 2 \ 2 + \w 3 \ 2 ?1 — ^2 l+\w 1 \ 2 +\w i -w 1 w 2 \ 2 



-i W2 — W3 — W2+W3 p —i (l+Wi)(w3-WiW2)-(l+Wl)(w3-WlW2) 

y/2 1+\W2\ 2 + \W3\ 2 ^ 2 ~ \/2 l + |«Ji| 2 + |?«3-?«iUJ2| 2 

\w 2 \ 2 -\w 3 \ 2 c _ l-|™i| 2 



^ 3 1 + |«)2| 2 + |«J3| 2 \+\w\\ 2 + \w3-W\W2\ 2 ' 

/■ • W 2 W3— W2W3 £ • W\— W\ 

S,4 — , |„,„|2 r |„,,|2 <~,4 — ' ill,.,. 12 , |,„,_,„„„J2 5 



' l+|«J2| 2 + |w>3| 2 ^ 4 l+|lUi| 2 + |lU3— w\w 2 \ 

/• _J_ W2 — W3+W2 — W3 p 

^5 ,/H 14-|, /10 |2 i |„,„|2 V 



1 W2-W3+W2-W3 £ 1 (1+W 1 )(U)3-W 1 W2) + (1+W 1 )(W3-W 1 U)2) 

y/2 1+|«>2| 2 + |«)3| 2 ?5 — ^ l+|l«l| 2 +|w)3-t«l«!2| 2 



i W2+W3—W2—W3 t - j (l — 1U 1 )(w3 — W 1 W2)—(l—Wi)(w3—WiW2) 

^ 6 ~ V2 1+M 2 +|">3| 2 ?6 ~~ V2 l+JwiP + l^s-^l^l 2 ' 

/• 102^3+^2^3 C W\+W\ 

S>7 — ij.|,,,„|2 . |,„ i2 



l+|l02| 2 + |«'3| 2 ^ 7 1 + |u)i| 2 + |ui 3 -«)iUI2| 2 ' 

> _ 1 2-|w 2 | 2 -|i"3| 2 p _ 1_ l+|wi| 2 -2|tt)3--wit»2| 2 

^ 8 l + M 2 + h 3 | 2 ^ 8 ^ l+hl| 2 +|w3-WlW 2 | 2 " 

Here Wi, are complex parameters on a generic orbit, m and g are initial values 

of //3 and us respectively. The initial values fix an orbit. For a degenerate orbit one has 
to assign m = and wi — 0. 

After this parameterization the effective Hamiltonians get the form 



fc=l a,p 



deg _ \^ ^Ca ^Ca 

dw a dW/3 ' 



deg _ V- 



gen 

91b 



\ Hlllr, Fill) FilDr, 



dw a dwp dw a dwp dw a dwp 

The tensors g gen and g deg serve as metrics on orbits in terms of the complex parameters 
w = {wi, wi, W2, W2, W3, W3} for a generic orbit, and w = {w 2 , w 2 , W3, W3} for a 
degenerate orbit. Note, that the metrics do not depend on the initial values m and q, 
fixing an orbit. All generic orbits have the same metrics, as well as degenerate orbits. 



4. Results and discussion 



Our main result is the following. For a magnetic system with the spin s^l we propose 
two effective classical models that describe fluctuations of the mean field by the Landau- 
Lifshitz like equations. We consider the 8-component mean field fi = {fi a }^ =l , taking 
into account not only magnetization but also quadrupole moments. 

The effective models deal with large-scale (slow) fluctuations of the mean field. 
Small-scale (quick) fluctuations are cut off by quasiaveraging. In this process some 
observed quantities become rigid and serve as constrains determining the manifold where 
the mean field lives. This manifold appears to be a coadjoint orbit of the group SU(3). 

Also we propose a complex parametrization for the manifold and reduce the mean 
field and the Hamiltonian density to complex parameters. Remarkably, that in terms 
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of the complex parameters the density becames independent on boundary values of fi. 
Moreover, the Hamiltonian density serve as a Riemannian metrics on the manifold. 

In the case of SU(3)-invariant model we deal with the magnet whose ferromagnetic 
and nematic states at zero temperature are equiprobable. That is why we propose two 
effective Hamiltonians: "H gen for states with the ferromagnetic order at zero temperature, 
and %deg for states with the nematic order (when magnetization is zero) at zero 
temperature. Also we produce equations f }2"2"]) and f l2H]) describing large-scale fluctuations 
of the mean field [i near a nematic and a ferromagnetic ordered states respectively. 

The proposed classical models can be used to construct topological excitations [11] , 
which are stationary solutions of the Landau-Lifshitz like equations. These excitations 
realize destruction of a long-range order in 2-dimensional spin systems at nonzero 
temperatures, according to the Mermin- Wagner theorem. 

The considered scheme is easily extended to the case with higher powers of exchange 
interaction. For an arbitrary spin s the spin operators {S%} are defined over the 
(2s+l)-dimensional space of representation of the group SU(2). The complete matrix 
algebra generated by the spin operators is Mat( 2s +i) X (2s+i)- Then one can consider a 
spin Hamiltonian with powers of exchange interaction up to 2s. Such Hamiltonian 
admits a bilinear form, if one takes into account multipole moments. In the mean field 
approximation this quantum model corresponds to a Hamiltonian system on a coadjoint 
orbit of the group SU(2s+l). Each orbit has a Hamiltonian system, which serves as an 
effective classical model. 
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